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We address the possibility that CP-violation in Bs — Bs mixing may help explain the origin of 
the cosmic baryon asymmetry. We propose a new baryogenesis mechanism - “Electroweak Beauty- 
genesis” — explicitly showing that these two CP-violating phenomena can be sourced by a common 
CP-phase. As an illustration, we work in the Two-Higgs-Doublet model. Because the relevant CP- 
phase is flavor off-diagonal, this mechanism is less severely constrained by null results of electric 
dipole moment searches than other scenarios. We show how measurements of flavor observables by 
the D0, CDF, and LHCb collaborations test this scenario. 


Introduction The baryon asymmetry of the Universe 
(BAU) has been precisely measured by the WMAP col- 
laboration. Combining its five year results with those 
from other CMB and large scale structure measurements 
gives Qh? = 0.02265 + 0.00059 [1] which is in excel- 
lent agreement with the 95% C.L. range 0.017 — 0.024 
obtained from Big Bang Nucleosynthesis [2]. The im- 
plied ratio of baryon density ng to entropy s is ng/s = 
(8.82 + 0.23) x 107-1. 


To generate the observed BAU, three Sakharov crite- 
ria [3] must be satisfied in the early Universe: (1) baryon 
number violation; (2) C and CP violation; (3) a depar- 
ture from thermal equilibrium (or CPT violation). These 
requirements are not unconquerable, though doing so re- 
quires physics beyond the Standard Model (SM) of parti- 
cle physics. Indeed, there exist a number of possibilities, 
though none have been conclusively established. One of 
the most popular — standard thermal leptogenesis — pro- 
vides a theoretically attractive solution, yet it is generally 
difficult to test experimentally. It is, therefore, worth 
considering scenarios that may be more directly tested 
laboratory experiments. 


A particularly interesting and largely unexplored pos- 
sibility involves CP-violation that enters the B, meson 
system. The relevant phases are generically flavor off- 
diagonal, making them less susceptible to constraints 
from searches for permanent electric dipole moments 
(EDMs) that challenge other baryogenesis scenarios (for 
an illustration in the Minimal Supersymmetric Standard 
Model (MSSM), see, e.g. [4]) . Moreover, recent mea- 
surements in B-factories and at the Tevatron exhibit in- 
dications of CP-violation that differ by a few standard 
deviations from the SM predictions [5, 6], though even 
more recent results from the LHCb collaboration favor 
smaller effects [7]. From a theoretical perspective if the 
CP phase(s) encoding the CP-violation in B, system can 
successfully drive the generation of the BAU and can be 
probed experimentally, our understanding of the BAU 


problem will be considerably advanced. 

In this letter, we report on an initial effort addressing 
this question. We propose a new mechanism in the frame- 
work of Electroweak Baryogenesis (EWBG), explicitly 
showing that the CP-violating phenomena characterized 
by different energy scales (B, observables and BAU) can 
be sourced by a common CP-phase. As an illustration, 
we will work in a Two-Higgs-Doublet Model (2HDM). 
In this context, if a sufficiently strong, first-order elec- 
troweak phase transition (EWPT) occurs in the early 
Universe, the CP phase associated with the tree-level, 
Higgs-b-s interaction at the phase boundary can induce 
CP-asymmetries that ultimately induce the BAU. While 
EWBG in the 2HDM has been discussed extensively in 
the past [8], including two recent studies using the 2HDM 
that have addressed the possible connection between the 
BAU and B, observables|9] (see also [10]), we empha- 
size that the mechanism discussed below is the only one 
thus far to explore the feasibility of baryogenesis directly 
driven by the b > s CP-violation. Given the novel fea- 
tures that are generically absent elsewhere and the crucial 
role played by “beauty” quarks, we denote this mecha- 
nism “Electroweak Beautygenesis” (EWBTG). 

In what follows, we concentrate on the issue of CP- 
violation and do not treat the question of the first order 
EWPT in the 2HDM. Following Ref. [9], we instead re- 
fer the reader to more general studies that may indicate 
its feasibility [8]. We note, however, that these anal- 
yses are typically gauge-dependent and therefore open 
to question. Rather than delve into these subtleties of 
perturbative treatments of the EWPT, we also refer the 
readers to a recent and more extensive discussion [11]. 
Two Higgs Doublet Model The Higgs sector in the 
general 2HDM is (Hy,4 are Higgs doublets with their SM 
gauge charges being (0, 2,+1/2)) 


L = MO (eH uh — AL. Q' Hadh 


-yQ Huuh + yi Q' (Hidh + h.c.. (1) 


In a supersymmetric embedding, the first term can arise 
at loop level [12]. For experimental relevance, we focus 
on the two flavor b-s system, with its mass matrix 


Mij = YijVu + AijVva = (Yaz sin B + Aij cos B)u, (2) 


where vu, a are vacuum expectation values (VEVs) of the 
neutral Higgs fields with v = \/v2+v% and tang = 
Uy/Ud- Vua are functions of spacetime during the EWPT. 
Meanwhile, Hy, = — cos GH, + sin BH}, a linear combi- 
nation of Higgs mass eigenstates, will introduce flavor- 
changing neutral current effects at zero temperature. 

Since we are investigating the feasibility that a com- 
mon phase can source the BAU and account for the Bs 
CP-violating observables, we will work in a simplified 
but sufficiently representative scenario, deferring a more 
comprehensive treatment to future work where the fol- 
lowing scenario would arise in one region of parameter 
space. First, we take tan 8 = 1 at zero temperature, em- 
phasizing that tan 8 is not a constant during an EWPT. 
Second, we assume Ysb = Asp = Msp = 0. In the limit 
of Yss, Ass — 0, there is one CP-violating phase after ap- 
propriate field redefinitions. Without loss of generality, 
we assume that Aps is complex (with 0),, = Arg(Abs)) 
and yos, Yob and App are real, and furthermore, assume 
Aii = Yii and |Ass| = |yos|. The mass matrix is then 


+2... 0 
Ge + etos) oe) a (3) 


here &; = |Ai;|/V2 and the “+” signs are due to Yss, Ybs 
and yy, > or < 0. Denoting mp, as Mas = Aexp(i@), we 
have A = 2&.| cos(@),, /2)|v, 9 = @,,/2 for yos > 0, and 
A = 2&).| sin(@y,, /2)|v, 6 = (Oar. + 7)/2 for Yos < 0. 

The mass matrix can be diagonalized by a unitary 
transformation Ut MUR = D. In the small m,, limit, 
U_ is trivial and Up is parametrized by a rotation angle 
ar = —arctan(A/mpp). The coupling of Hss and br, sR 
quarks in the mass eigenstate basis is given by 


dos = S0s[-F1 + exp(iy,, )] cos aR. (4) 


with Arg(¢),) = @+7/2 for Yes > 0 and < 0, respectively. 
It is just the phase 0 (or 0),,) that both introduces CP- 
violation in b — s transitions (via Cs) and source the 
generation of baryon asymmetry (via Mas). 

Electroweak Beautygenesis Production of the BAU 
during a first-order EWPT involves a dynamic generation 
of CP-violating charge asymmetries through particle in- 
teractions in the wall of nucleated bubbles. Those charge 
asymmetries are converted, via left-handed fermions 
(nz), into the baryon asymmetry through the elec- 
troweak sphaleron transitions. We ignore the wall cur- 
vature in our analysis so all relevant functions depend 
on the variable Z = z+ vwt. Here vy, is the wall ve- 
locity; Z < 0 and > 0 correspond to the unbroken and 
broken phases, respectively; and the boundary extends 


over 0 < Z < Ly. As pointed out in [19], the trans- 
port properties of particles during the EWPT is most ap- 
propriately treated using non-equilibrium quantum field 
theory. Working in its closed time path formulation (for 
pedagogical discussions, see [20]) and under the “VEV- 
insertion” approximation (see, e.g., Refs. [9, 19-21]), we 
compute the CP-violating source induced by the Higgs 
mediated process bz —> sr — bz. It is given by 
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Here np(x)=1/(exp(x) + 1) is the Fermi distribution; 
Ebt sR = Wb, sr — Ilbz,sp are complex poles of the spec- 
tral function with wĝ, sp = k? + Mẹ, spi and Me, sp and 
Ty,,s2 are thermal parameters. This source corresponds 
to the “A”-type terms in Eq. (58) of [20], after properly 
accounting for temperature-independent vacuum contri- 
butions that are removed via normal ordering [31]. The 
quantity 6(2) = d0(z)/dt is given by 


io = A sienen 


with f (Z) = (vu(Z)va(Z) — vu (Z)ùa(Z)) ~ Vwv?6B/Ly be- 
ing a function describing the relative variation of the 
Higgs VEVs across the bubble wall. Although analy- 
ses performed in the MSSM [22] indicate 58 ~ O(107?), 
a systematic analysis is absent in the 2HDM. Here, we 
will simply adopt 68 = —0.05 (if ybs is complex and Xp, is 
real, we need 6{ > 0 to keep the sign of f(z) unchanged). 
Note, S id is non-zero only within the moving bubble wall, 
where 6(Z) Æ 0. 

In contrast to EWBG driven by flavor-diagonal 
sources, the transport of both the second and third fam- 
ily particles is sourced by CP-violating terms. We de- 
fine the number densities {Q1,23,U,D,C,S,T,B,H = 
Ht + H} — H} — HÌ} which correspond, respectively, 
to left-chiral quarks of different families, right-chiral up, 
down, charm, strange, top and bottoms, and Higgs 
bosons. Since all light quarks (except bz and sp) are 
mainly produced by strong sphaleron processes and all 
quarks have similar diffusion constants, baryon num- 
ber conservation on time-scales shorter than the inverse 
electroweak sphaleron rate implies the approximate con- 
straints Qı Qo 2U 2D 2C 2B and 
S +T +Q = 0. The set of Boltzmann equations is 


I Q3, = Im, (Er = Soy) +T: (Er = EH — Eqs) 
+20 ss (Er — 2g, + £s + 8EB) + SF 
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Here 04 = vw $ — Data in the planar bubble wall 
approximation with D, being a diffusion constant, and 
Ea = Na/ka with na and ka being the number density 
and the statistical factor of particle “a”. Apart from the 
CP-violating sources, the interactions in Eq. (7) include 
(i) inelastic top Yukawa (T,) and strong sphaleron (Tss) 
processes; (ii) top relaxation processes (Tm,), while we 
neglect the other Yukawa interactions since T D4/v2, < 1; 
and (iii) Higgs relaxation processes (T) due to Higgs 
mass mixing, with typically T'a < I'm, in this scenario. 
Assuming SZ (z < 0) = 0, we solve the Boltz- 
mann equations for the net left-handed fermion density 
nr, = y Qi analytically order-by-order in 1/T ss, with 
Tss = 16047. The leading contribution arises at first or- 
der in this expansion. The baryon asymmetry pp is then 
produced in weak sphaleron process, described by [23] 


15 
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where Tus = 120a3,T is the weak sphaleron rate [24]. In 
the broken phase this gives (ks = kp is assumed) 
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Vw)/2D. Here D, T and S are, respectively, the effective 


diffusion constant, decay rate and CP-violating source 
for the Higgs number density. D is defined in [20], while 


with r = 


T = (9kr +9kg +kp)(Um, + 20n)/X 
5 = kulkos — Tkr + kp) SOX (10) 
KX = 9ka,kr + kBkas +4krkpt+ ky (9kr + 9kas + kp). 


Note that while the weak sphaleron transitions are driven 
by the diffusion tail for nz that extends ahead of the ad- 
vancing wall in the unbroken phase (Z < 0), the solution 
in Eq. (9) contains an integral over the source in the bro- 
ken phase that appears when matching the solutions to 
the Boltzmann equations at the phase boundary. 
CP-violation in the B, — B, Mixing Depending on 
the details of the scalar potential, Hy, may be approxi- 
mately a mass eigenstate, which we assume for illustra- 
tion. Tree-level exchange of Hy, with a VEV insertion 
leads to a B, — B, mixing operator in the basis of quark 
mass eigenstates (Aj; is an effective new physics scale): 
Cia 


re (bLsR)(bLSR), with mia, ~ UNbs 


bs 


The RG running of this operator involves a mixture of 
scalar operator Ofer = (bPrs)(bPrs) and tensor op- 
erators Of'gR = (bo"” Prs)(bo pu Prs) whose matrix el- 
ements are (B,|O%,p|Bs) ~% —5mp, f2, Bogp/24 and 
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FIG. 1: Contours of ng/s (in units of 107'°). Orange, yel- 
low and green contours indicate 95% C. L. Tevatron (with 
and without di-muon asymmetry) and LHCb constraints, re- 
spectively. We assume kg, = 2kr = 2kg = 6 [25], ky = 4, 
Vw = 0.4 [26], Ly = 2/T, D = 6/T and Ay, = 1 TeV. 


(Bs |\Oe RRBs) © mp, fh, Be pp/2 [14]. Assuming m$, 
~ (m, + m)? and BY, ~ Bier = Bp,, we ob- 
tain Mf, = (B.|H|Bs) = —G,f2,me, Ba, (5nsrr/24+ 
nrrR/2)/Ap,, with NSRR~ 1.87, 1TRR~ —0.01 [14]. 


gm sail fa,./ Be. 
0.019 + 0.001 0.27 + 0.15 (275 + 13) MeV 
(Ame, )™ (Ams,)*” [17] oe” 
19.304+2.2 ps7} 17.7740.12 ps7! (4.2 1.4) x 107° 
(Aly (ATs) 7° [5] A? [6] 
0.098 + 0.024 ps~* 0.097 + 0.031 ps~! (—7.40 + 1.93) x 107° 
(Amp, ) "SP [7] (ATR? [7] Bice [7] 
17.725 + 0.049 ps~! 0.123 + 0.030 ps4 —0.015 + 0.087 


TABLE I: The theoretical input parameters [2, 15] and the 
experimental data from the Tevatron and LHCb. 
Choosing T4 to be real and parametrizing Mf, as [16] 


M3, = (M3,)™A, with A, = |A,le’* , we have 
AT, = ATM cose +05); Am, = Am“ |A, 
_ ATM sin(ge™ + 2) 
“SL Am™M As. 
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26, = 28)" =F (11) 


Here Am, and AT, are the mass and decay width differ- 
ence between the heavy and light B, mass eigenstates, 
ag, is the charge asymmetry in semileptonic B, decays, 
and 8, measures the time-dependent CP asymmetries in 
the hadronic B, decay. 

The theoretical inputs and experimental results are 
listed in the Table I. The decay constants and bag pa- 
rameters are taken from Ref. [18], while A2,, and BPP, 
AT? are obtained by combining the DØ and CDF mea- 
surements [5, 6]. We perform a x? fit to the four ob- 
servables in Eq.(11), neglecting the correlation between 


AT, and $8, for simplicity. Assuming A,, of 1 TeV, we 
scan over the remaining parameters, yielding the regions 
of 95% C. L. from the Tevatron and the LHCb results. 
Discussion The contours of constant ng/s in the 
sign(Yyos)&bs — Sin@,,, plane are indicated in Fig. 1. We 
observe that the regions favored by the low-energy fla- 
vor studies at 95% C.L. overlap with regions in which 
a sizable portion of the baryon asymmetry is generated. 
The LHCb results on the B, hadronic decay are more 
constraining on the parameter space than the Tevatron 
ones that do not include the dimuon asymmetry. Al- 
though tension exists between the LHCb- and Tevatron- 
favored regions, it appears feasible that a common CP- 
violating phase may be responsible for both generating 
a non-negligible portion of the BAU and accounting for 
observations in the B, system. 

A definitive statement awaits the resolution of both 
the experimental tensions as well as several theoretical 
issues, including the development of a VEV-resummed 
CP-violating source (for recent progress, see, e.g., [27]), 
analysis of the full numerical solutions to Eqs. (7), and 
completion of a gauge-invariant analysis of the EWPT in 
the 2HDM. Indeed, the results of this initial study are 
likely to indicate the maximum magnitude of the BAU 
that can be achieved in this scenario, given the generous 
assumptions we have made about various input param- 
eters, including 66 and vw and the use of an analytic 
rather than numerical solution of the Boltzmann equa- 
tions. Nevertheless, we expect that after future refine- 
ments are implemented, EWBTG may account for an 
interesting portion of the BAU in appropriate regions of 
parameter space. It thus appears that a more compre- 
hensive analysis of this scenario is warranted. 

Though the foregoing discussion relied on the illustra- 
tive case of a two-flavor system of the 2HDM with a sin- 
gle phase, generalization to variants, including minimal 
flavor violation with flavor-blind phases (e.g., see [28]) 
or spontaneous CP-violation, would be straightforward. 
We leave the consideration of these possibilities, along 
with EWBTG in other models such as the four-family 
SM (e.g., see [29]), family non-universal U(1)’ model [30], 
and supersymmetric models, etc. to future work. 
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